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ABSTRACT 


In applying the charge exchange mechanism to ion phenomena within 
the Earth's magnetosphere it is critical to the proper interpretation 
of observations that the charge exchange lifetimes for the ions be known 
as accurately as possible. Various new results have been published which 
significantly modify the charge exchange lifetimes which have been used 
in space physics research during the past decade and a half. Some of 
the newer results have been used in the application of the charge exchange 
decay mechanism but the use has been limited and for the most part incomplete. 
The neutral hydrogen density distribution now yields lifetimes which are 
shorter than previously calculated, while the functional dependence of the 
lifetimes on pitch angle provides for slower decay for ions mirroring off 
the geomagnetic equator. This review coalesces and summarizes the latest 
and best measurements of the physical quantities involved in the complete 
calculation of the charge exchange lifetime of the mirroring magnetospheric 
ions. 


INTRODUCTION 


The exosphere of atomic hydrogen surrounding the earth extends far 
deeper into the magnetosphere than it was believed in the early sixties. 

The density of atomic hydrogen has not yet been measured directly in the 
magnetosphere, unfortunately, and an uncertainty in our understanding of 
associated geophysical phenomena therefore remains with us. Fortunately, 
however, the neutral hydrogen density can be inferred from Lyman-alpha 
emission measurements (Bertaux and Blamont , 1973; Vidal -Mad jar and Bertaux , 
1972; Meier and Mange , 1970) or from ion mass spectrometer measurements 
(Brinton , 1975), coupled with the theoretical models for radial distance 
dependence (Johnson and Fish , 1960; Chamberlain , 1963). The derived 
neutral hydrogen density is large enough so that the process of charge 
exchange in the trapping region, where the energetic protons and heavier 
ions can be considered lost through the acquistion of one or more charge- 
neutralizing electrons from the ambient hydrogen, is a phenomenon to be 
seriously and exactly considered. Indeed, since the early suggestion of 
Dessler and Parker (1959) that the charge exchange decay mechanism could 
be significant in the ring current decay, arguments for and against this 
mechanism have been brought forth by various authors (Frank , 1967; 

Swisher and Frank , 1968; Prblss , 1973; Smith et al . , 1976; Tinsley , 

1976; Lyons and Evans , 1976). Smith et al . , (1977) present several 
aspects of these arguments and draw conclusions about the ionic composition 

O 

of the ring current based on Explorer 45 (S -A) observations. 

In virsw of the importance of charge exchange decay as a loss mech- 
anism for magnetospheric ions, this paper coalesces and summarizes what 
are the latest and best measurements of the physical quantities involved 
in the calculation of the charge exchange lifetime of the mirroring ions. 
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These quantities are the neutral atomic hydrogen density, the charge 
exchange cross section of various ions and the lifetime dependence of 
mirroring ions on the mirror latitude. 

There are two principal relations for the charge exchange mechanism 
in the magnetosphere: 

^e ~ n(r )av ^ 

0 


where, is the mean lifetime of protons or any other ion species con- 
fined to the equatorial plane for charge exchange decay with atomic hydro- 
gen, n(r^) is the hydrogen density in the equatorial plane, v is the 
velocity of the ions and a is the charge exchange cross section of the 
ion with atomic hydrogen, and 


m 




cos'^A 


m 


( 2 ) 


where, is the charge exchange lifetime of ions mirroring off the geo- 
magnetic equator at the mirror latitude. 

The early work by Liemohn (1961) used a and n(rj^) measurements avail- 
able at that time to provide quite complete graphs for as a function 
of proton energy and the radial distance of the protons. He also obtained 
T as a function of mirror latitude in the approximate form of x^ = x. cos^A 
Until recently, these results have been the basis for all discussions of 
the charge exchange mechanism in the magnetosphere, even though better 
measurements and theories have existed for both n(r^) and a. Smith et al . 
(^'976) were one of the first to use charge exchange lifetime values modi- 
fied from those given by Liemohn (1961) in explaining satellite observations 
by using better determinations of n and a in equation (1). Tinsi ey (1976) 
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also modified x_ by usinq updated values of n{r„) and a and Lyons and 
0 0 

Evans (1976) used these values in explaining their observations. Smi th 
and Bewtra (1976) made an additional step forward by recomputing the 
relationship between lifetimes of particles mirroring off-equator to the 
lifetime of particles confined to the equator (equation 2). They deter- 
mined that the relationship was 


_ ^ 3.5+0.2 

T = T cos \ 

me m 


and, therefore, of f-90° pitch angles fluxes would not charge exchange 
decay as rapidly as was previously thought. Independently Cowl ey (1977) 
in doing a similar calculation art'ived at a nearly identical conclusion. 


GEOCORONAL HYDROGEN DENSITY MODELS 


The neutral density distribution in a planetary atmosphere is a 
function of the boundary conditions at the exobase (i.e., the tempera- 
ture and density distributions at the exobase). The exosphere is 
defined to be the region of a planetary atmosphere where the number of 
atomic collisions per unit time is negligible and where the controlling 
factors are gravitational attraction and thermal energy conducted from 
below. Thus below the exobase (also called the critical level) it is 
generally assumed that the constituents have an isotropic Maxwell- 
Boltzmann distribution, and above this level, in the collisionless 
medium, one uses Liouville's theorem to get the exospheric density as a 
triple integral of the distribution function over the momentum space. 

In this exosphere particles with an upward velocity will describe 
ballistic trajectories. Depending on their velocities at the exobase, 
they will either describe hyperbolic trajectories and escape from the 
atmosphere, or elliptic trajectories and return to another point on the 
exobase. 

Johnson and Fish (I960) formulated a model for hydrogen density 
which had been extensively used in the sixties and early seventies for 
the study of magnetospheric problems. They hypothesized that in the 
exosphere the particle distribution corresponds to a complete Maxwellian 
distribution, except with certain types of hyperbolic trajectories ex- 
cluded. They also assumed the conservation of escaping flux, and obtained 
a numerical solution to the hydrogen distribution in the exosphere. In- 

n 

dependently, Opik and Singer (1959, 1960, 1961) gave an analytical 
expression for both types of particles, elliptic and hyperbolic. Cham- 
berlain (1963) added a third class of particles, namely, satellite 
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particles, created by a small number of collisions in the exosphere. In 
Chamberlain's terminology, these particles correspond to elliptic orbits 
of particles having perigee outside the critical level. Thus, these 
particles do not correspond to trajectories originating from within the 
exobase. On the other hand, the particles with elliptic orbits having 
perigee below the critical level rise from the critical level and even- 
tually fall back. These are called ballistic particles. The third class 
of particles correspond to hyperbolic orbits and are called escaping par- 
ticles. Chamberlain also considered the sources and sinks of geocoronal 
particles- He examined the equilibrium questions that arise because of 
photoionization or other loss mechanisms and related this to fixing the 
abundance of satellite particles. This was done by introducing the con- 
cept of a satellite critical altitude Below this altitude there 

exists a complete isotropic distribution of energetically allowed satellite 
particles, and there are no satellite particles which have a perigee above 
this altitude, due to the assumption that the collisional creation of 
satellite orbits is negligible. 

The models described above correspond to an exospheric density 
distribution with spherical symmetry i.e., atmospheric structure is 
assumed uniform over the globe and there are no latitude or day-night 
effects. However, the interpretation of ground based, rocket and satellite 
experimental results often requires some degree of asymmetry in the 
exospheric density distribution of hydrogen ( Meier 1969; Metzger and 
Clark , 1970). Vidal -Madgar and Bertaux (1972), McAfee (1967) and 
Quessette (1972) have extended Chamberlain's arguments and introduced 
axial asymmetries. Vidal-Madjar and Bertaux have an involved model 
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requiring numerical integration corresponding to any temperature and 
density distribution at the exobase. They derive density at a particular 
point in the exosphere by integrating the distribution function over the 
momentum part of the phase space. However, all the following discussions 
shall be confined to spherically symmetrical models only. 

The density n{r) for the Chamberlain model is expressed in terms of 
density at the exobase, N^, and a partition function, as follows: 

n(r) = i;(X) 

where X(r) = 

G = Gravitational constant 

= Planetary mass 

M = Atomic mass 

k = Boltzman constant 

T^, = Temperatur'^ at the exobase 

For neutral hydrogen, X^can be expressed by assuming the exobase to be 
at 500km by 

, _ 6981.4 

^c " T^ (Tn K) 

and 

wpi . 1-0785 Xp 

^ r (in Earth Radii) 
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The partition function, c{X) , can be calculated for the appropriate 
type of particles by integrating the distribution function with proper 
boundary conditions of moments. We summarize the partition functions 
for the ballistic, satellite and escaping particles as: 

*^bal ^ 


r (X) = e-'f'l Y (3/2,X-,j) ) 

sac ^ * 


W (X)-^[#-y(3/2.X) 


- e"*l - T (3/2,X-t|.,))j 


where 


’1^1 = 


1 X + X, 


and Y is the incomplete r- 


f unction (Abramowitz and Stegun , 1965), 

As was stated earlier, Chamberlain introduced the concept of satel- 
lite critical altitude, which will assume a role for satellite par- 
ticles analogous to that of cricital level, R^, for ballistic and escaping 
particles. Thus, for distances below R^^ there is, because of occasional 
collisions, a complete isotropic distribution of the satellite particles 
and the partition function for r < R^^, is 

^ ~ ^bal ^ ^sat ^ *”esc 

Y (3/2.X) + (X) 

For, r > Rg^, only those satellite particles are present which have 
perigees between R^ and Thus, the ballistic formula can be modified 

as though these orbits arose from R rather than R . (The ballistic or 
escaping components are not actually affected by this modification). 


Thus for r > R^^, the partition function is 

^ " ^esc ^bal ^^»^cs^ 
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( 5 ) 


We see that corresponds to using only equation (5) and 

implies no population of satellite particles. On the other hand, R^^= «■ 
corresponds to using only equation (4) and implies a complete distribution 
of satellite particles. A realistic situation, however, corresponds 
to a mixture of these cases and most Lyman-alpha measurements are ade- 
quately described by a value of R = 2.5R ( Carruthers et al . , 1976; 

w ^ L» 

Meier and Mange , 1970). 

In Figure 1, we show the functional dependence of atomic hydrogen 
density on the radial distance from the center of earth as given by the 
Johnson and Fish (1960) model (solid line) and Chamberlain (1960) model 
(dashed line). A unit density at ^-he exobase (500 km) and the value of 
T = 1200°K are assumed. Additionally, for the Chamberlain model R = 

2.5 R^ is assumed. As is obvious, the Chamberlain model gives much 
higher densities than the Johnson and Fish model. In Figure 2, we show 
the atomic hydrogen density distribution as given by the Chamberlain 
model for four exobase temperatures, T^, (900°K, 1000°K, 1100°K and 
1200°K). For a specific exobase temperature, a unit density at the exo- 
base is assumed and the density distributions for 1) no satellite particle 
distribution, 2) partial population of satellite particles (correspond- 
ing to Rg^ =2.5 R^) and 3) complete satellite particle distribution are 
shown. 

Thus, for a known atomic hydrogen density and temperature at the exobase 
the density at an observation point can be easily determined from these 
graphs. 

In the absence of a simultaneous measurement of density at the exobasc 
and the magnetospheric observation, one has to employ empirical models 
to derive the exobase paiameters required by „:ie Chamberlain model. 
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Hedin et a1 . , (1977) have developed a global thermospheric model (MSIS 


model) for the neutral thermospheric composition and neutral temperature 
based on mass spectrometer data from five satellites (AE-B, OGO- 6 , San 
Marco-3, AEROS-A, AE-C) and incoherent scatter measurements at four 
ground stations (Arecibo, Jicamarca, Millstone Hill, St. Santin). The 
input parameters for this model are the altitude, latitude, longitude- 
seasonal parameters, magnetic activity parameter (Ap index), and the solar 
parameters (daily and average 7 flux). It should be noted that the 
temperature and density distribution over the exobase can be represented 
by an approximate distribution having an axial symmetry and nearly sinu- 
soidal variation between minimum and maximum. Thus, the average exobase 
temperature can be effectively approximated by J 5 

worth pointing out that the hydrogen density distribution is inverse in 
phase with the temperature. Thus, the actual variation of the atomic 
hydrogen (see Figure 2) in the magnetosphere is greatly reduced. Vidal - 
Madjar et al . (1974) have parameterized the atomic hydrogen density and 
temperature at the exobase based on the OSO- 5 solar Lyman alpha measurements. 
They provide in a tabular form for the years 1969 to 1971, the daily 
exobase temperature (including minimum and maximum) and the atomic 
hydrogen density as well as the fit through their data showing relationship 
between exobase temperature and the density. In Figure 3, we show the 
actual range in the atomic hydrogen density variation in the magnetosphere 

for a variation in exobase temperature of 900°K to 1200°K and the cor- 

. . . ■ .4-3 5-3 

responding variation in the density of 4 x 10 cm to 1.1 x 10 cm . 
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CHARGE EXCHANGE CROSSSECTION MEASUREMENTS 


As was stated earlier, energetic protons and heavier ions can be 
lost from the trapping regions through acquisition of charge-neutralizing 
electrons from the ambient neutral hydrogen surrounding the earth. This 
is the charge exchange process. Charge exchange lifetimes are sufficiently 
short at the ring current energies that this process is an important 
loss mechanism for these ions. In this section, we shall summarize the 
survey of cross sections, measured and theoretically derived, for the 
following collision processes: 


H"^ + H 

— 

H + H’*' 

(a) 

+ H 

— 

*4“ 

He + H 

(b) 

He’*' + H 


He + H’*' 

(c) 

O'*" + H 


+ 

0 + H 

(d) 


In these interactions, it is assumed that each member of the col- 
liding pair is in its ground state and the final products may be in the 
ground state or in an excited state. Momentum and energy transfer 
during charge exchange are negligible, i.e., the incident particle 
leaves the collision with essentially unchanged energy but a changed 
charge state. This is true whenever the incident energy is large compared 
with excitation and ionization energies, which are of the order of a few 
eV. It is also important to point out that ior many charge exchange 
cross section experiments, molecular hydrogen was used rather than 
atomic hydrogen as the target. In all the following discussions, the 
cross sections have been properly normalized to atomic hydrogen as the 
target. 

H~^ on H 

Cross sections for reaction (a) at low energies (less than 50 keV) 
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were first measured by Fite et a1 . (1960) and used extensively in early 
work on the studies of ring current particle loss. The best experimental 
measurements available for the low energy region are those of HcClure , 
(1976) in the energy range of 2-117 keV. He used a heated target chamber 
containing mostly dissociated hydrogen gas. His estimat'dd error is ±5^ 
after correcting for the presence of molecular hydrogen in the target gas. 
Several theoretical calculations ( Mapelton , 1962; Bates and Dalqarno , 1953; 
Dalqarno and Yadav , 1953) agree with his measurements over different parts 
of the energy range. Fite et a1 . (1960) used a crossed beam experiment 
and their measurements should provide the cross section reliably in the 
low energy end (less than 10 keV). Stier and Barnett (1956) provided the 
experimental data for the energy range 120 to 200 keV using a molecular 
hydrogen target and give an estimated error of tS%. For completeness » 
the best data for high energies were provided by Toburen et al. (1968) 
in the energy range 200 keV to 1 MeV with an error estimate of ±10%. 

These values are summarized in Table 1. 

He^~*~ on H 

Fite et al . (1962) provided the experimental results for reaction 
(b) in the energy range 100 eV to 32 keV using an atomic hydrogen target. 

3 

They used modulated crossed beams and the isotope He in order to dif- 

++ + 

ferentiate the charge to mass ratio for He and H 2 . At the low energy, 

data are good to within +25% with the error decreasing to about +13% at 

high energies. They normalized their data to the charge exchange cross 

sections for protons on atomic hydrogen at each energy. Shah and Gil body 

(1974) have provided the best measurements for this reaction in the energy 

3 ++ 

range 6 to 60 keV by the passage of He ions through a tungsten tube 
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+ + 

TABLE 1. Charge Exchange Cross Section for H + H ^ H + H . 
See text for source of measurement. 


2 

Energy (keV) Cross Section (cm ) 


1.92 

13.3 

X 

3.04 

12.1 

X 10"^^ 

3.82 

11.1 

X 10“^^ 

4.8 

10.5 

X 10"^^ 

6.05 

9.85 

X 10“^^ 

9.6 

8.6 

X 10'^^ 

12.1 

7.5 

X 10"^^ 

T 9. 2 

5.0 

X 10-'^ 

30.4 

2.97 

X 10-'“ 

38.2 

1.86 

X 10-'“ 

48.0 

1,10 

X 10-'“ 

60.5 

6.0 

X 

o 

1 

76.2 

3.0 

X 10-'^ 

100.0 

1.24 

X 10-'^ 

117.5 

6.4 

X 10-'® 

120.0 

6.6 

X 10-'® 

140.0 

3.7 

X 10-'® 

160.0 

3.7 

X 10-'® 

180.0 

1.4 

X 10-'® 

200.0 

8.8 

X 10-'® 
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^ -J. 

furnace containing highly dissociated hydrogen and observed He ions in 
all final capture states. They estimate an uncertainty of 8 to 1055. 

In general, the cross sections obtained by Fite et al . (1962) are much 
lower than those of Shah and Gil body (1974). .The Fite et al . (1962) 
data for reaction (b) was renormalized by Shah and Gil body (1974) to the 
cross section for reaction (a) as measured by McCI ure (1966). These 
values then show reasonable agreement' to the similar cross sections 
measured by Shah and Gil body (1974). 

Several theoretical calculations have been done for reaction (b). 

Basu et al . (1967) and Malaviya (1969) obtained n=2 capture cross sections 
by considering four final states of He'*'. Rapp (1973, 1974) has performed 
the most complete calculation for capture and excitation in reaction (b) 
and has considered 11 final states. These calculations show good agree- 
ment with earlier mentioned measurements, particularly in the 10 to 60 keV 

^ »[«i 

incident He energy region. 

It may be worth noting that most experimental cross sections are based 
on the mass 3 isotope of He , while most theoretical cross sections are 
for the more abundant isotope corresponding to mass 4. Assuming the scaling 
of the cross section, that is, the cross section depending on the velocity 
of the incident ions, it is a linear shift in the energy scale of the cross 
section for one isotope compared to the cross section for the other isotope. 
We shall be discussing the isotope corresponding to mass 4. Table 2 gives 
cross sections for reaction (b) as determined from a smooth line which 
we constructed through the data values obtained by Shah and Gilbody (1974) 
and the renormalized Fite et al . (1962) values. 
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TABLE 2. 

4 ++ 

Charge Exchange Cross Section for He 
See text for source of ineasurements. 

+ 

H V ^He^ + H*^. 

Ene. •' 

(keV) 

2 

Cross Section (cm ) 

5.5 


3.5 

X 

10’'^’ 

8.0 


5.5 

X 

lo-'^ 

12.0 


7.2 

X 

10-'' 

15.0 


8.4 

X 

10-1^ 

26.0 


12.0 

X 

10-1*' 

40.0 


12.0 

X 

lo'l*' 

52.0 


11.2 

X 

10"'*’ 

65.0 


10.5 

X 

lO-l*’ 

80.0 


9.0 

X 

10"'*’ 
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He^ on H 

All 1 son (1958) and Barnett and Stier (1958) have provided the measure- 
ments for this cross section at low energies and high energies, respect- 
ively. Allison and Garcia-Munoz (1962) have provided a report on the 
results of Allison (1958) and Barnett and Stier (1958). All measure- 
ments of this cross section have been with molecular hydrogen, so that 
results at lower energies will be relatively inaccurate. A possibility 
exists that the cross section for reaction (c) can be obtained by applying 
the principle of detailed balance on the inverse reaction; H + He H + He . 
However, the capture of He^ into an excited state for reaction (c) can be 
comparable to the capture into the ground state in certain energy ranges. 
Consequently, the inferred cross section may not be accurate. In Table 3 
the cross sections for reaction (c) which we obtained by a best fit of the 
values of Allison (1958) and Barnett and Stier (1958) are listed. 

o'*' on H 

Fite et a1 . (1962) have made the only measurement of the cross 
section for reaction (d) using a modulated crossed beam technique. 

However, in their experiment, the energy range of the incident 0 ions 
was 30 eV to 10 keV. Above this energy they have provided an empirical 
formula in a form which is characteristic of resonance behavior; the 
reaction (d) falls in this category due to the approximate energy 
balance. Fi te et a1 . (1962) use the cross section of reaction (a) at 
the incident ion energy of 1.9 keV to normalize the experimental data 
for reaction (d). As stated earlier, better measurements ( McClure , 

1966) for reaction (a) exist now, and this requires that the measurements 
of Fite et al. (1962) and their empirical formula be properly renormalized. 
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TABLE 3. Charge Exchange Cross Sect ion for ^He^ + H '^He + 
See text for source of measurement. 


Energy (keV) Cross Section (cm^) 


4.0 

3.6 

X 10 ‘^® 

8.0 

4.7 

X 10-18 

14.0 

6.7 

X 10"18 

20.0 

8.2 

X 10"18 

26.0 

9.6 

X 10-18 


11.6 

X 10"18 

50.0 

12.3 

X 10-18 

60.0 

13..0 

X 10-18 

70.0 

13.0 

X 

j 

o 

1 

CO 

80.0 

12.0 

CO 

1 

o 

X 

100.0 

10.8 

X 10"18 

120.0 

9.2 

X 10-18 

140.0 

7.7 

X 10"18 

160.0 

6.4 

X 10-18 

180.0 

5.6 

X 10"18 

200.0 

4.9 

X io ; i 8 


f 

t 


i ' 

i 

5 

t ' 
! 


i 
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We renormalized these values and* they are shown in Table 4 for incident energies 
up to 25 keV. At higher ion energies, however, the predicted cross section 
becomes unreliable. A better estimate of the cross section at these energies 
can be obtained by using the principle of detailed balancing on the reverse 
reaction. 

+ 0 H + o'*" (e) 

If we define by the cross section for the direct reaction and by 

the cross section of the reverse reaction (at same center of mass velocity), 

it can be shown ( Mapelton , 1972) that 

(2S^ + 1)(2L^ + l)o^ = 2(2S. + 1)(2L. + 1)0^ 

where and are the initial and final orbital angular momentum and 
and are the initial and final spin momentum. If we assume that 
the atomic oxygen, which is in the ground state ^P, ionizes to the ground 
state of 0"^ ions (namely, ‘‘S state), it can be easily seen that 

0^(O'^(-S) + H->0(3p) + H”^) = 9/8a^(H'^ + 0(^P)'^H + O'^'C’S)). 

Stebbinqs et al. (1964) have measured the cross section for reaction (e). 

Stier and Barnett (1956) have also measured this cross section, using 
molecular oxygen, with lower uncertainties. We have used these measure- 
ments and the principle of detailed balance to obtain the cross section for 
reaction (d) at 0 energies greater than 25 keV, and these are shown in 
Table 4. 

In Figure 4, we summarize the results for reactions (a), (b), (c), 
and (d) by plotting the 'normalized lifetime' i.e., ^ for various ions 
in the energy range 1 to 200 keV. 
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TABLE 4 


Charge Exchange Cross Section for the Reaction 0^ + H 0 + h’’^. 
See text for source of measurement 


Energy (keV) 


2 

Cross Section, (cm ) 



DEPENDENCE OF CHARGE EXCHANGE LIFETIME ON MIRROR LATITUDE 


The dependence of the charge exchange lifetimes on the mirror lati- 
tude had been considered by Lieniohn (1961) based on Johnson and Fish 
(1960) hydrogen density models. He obtained an empirical fit in the 
form of equation (2) with j=6. Since then, most of the researchers have 
been using these values to compare with their observations. Prdlss (1973) 
used a better hydrogen density model and speculated that j=5 would yield 
a better approximation. Smith and Bewtra (1976) and independently Cowl ey 
(1977) examined in detail the lifetime of mirroring particles and recomputed 
the relationship using the Chamberlain (1963) model of neutral hydrogen 
density. 

Smith and Bewtra (1976) and Cowley (1977) both considered that ions 
bouncing between the mirror points make a large number of bounces before 
charge exchanging with neutral hydrogen. Thus the instantaneous hydrogen 
density over a mirror path can be effectively averaged to get the time 
averaged density, n, that the particle encounters as 


n 



1 

^r(r) 


1 


ds 


ds' 


( 6 ) 


Using equation (1), then, t /t = n /n and the ratio /t can be computed 

> M C. J 1 1 

by calculating the integrals in equation (6). Smith and Bewtra (1976) 
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determined that, for = 1000°K, = 2.5R^ and at L = 4, a least squares 

fit of the calculated ratio, t„/t , to a function of the form cos'^X 

me m 

yielded a value of j = 3.4. Cowl ey (1977) also determined j to be in 
the range 3 to 4. 

The ratio is independent of N^, the density at the exobase, 

but does depend on the two other parameters, and R^^, of the Chamberlain 
model. Smith and Bewtra (1976) computed the effect on this ratio of vary- 
ing the two parameters. In fitting the ratio using various extreme combina- 
tions they found, for between 900°K and 1200°K, for R between 2.25R 

c sc c 

and 2.75R^, and for L-values between 2.5 and 10, that j fell in the range 
j = 3.5 ± 0.2. 

Cowl ey (1977) additionally suggested the use of an alternative 
approximation of the form 



(7) 


where A is the latitude at which the field line (cos^ A = 1/L) intersects 
the Earth's surface. He provides tabulations of Y3 for various L-values, 

T values and R values. He concluded that the alternative approximation 
(eq. 7) gave good values of over a large range and even better 

approximations than cos'^X^ and other functional forms which he tried. 

It is worth pointing out that at a mirror latitude of 35° the charge 
exchange lifetime using j = 3.5 is a factor of 2 longer than predicted 
by the Liemohn (1961) results. Another way of looking at this difference 
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is that using the new lifetime for protons of energy 10 keV mirroring 
at a latitude of 35°, the flux after 1 day will be a factor of 15 higher 
than that predicted by Liemohn's results. 
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SUMMARY 


In this paper we started with the two principal relations for charge .• 
exchange lifetimes. Figure 2 summarizes the normalized atomic hydrogen 
distribution as a function of radial distance using the Chamberlain model 
for a range of exobase temperatures and for various combinations of 
satellite particles. Thus, if the atomic hydrogen density at the exobase 
is known, one can determine the density at any point of observation. In 
general, the temperature and density at the exobase are functions of 
several geophysical parameters. Therefore, in Figure 3 we have shown 
the range in the actual density distribution for several magnetic storms 
during the early life of Explorer 45 in the period December, 1971 to 
May, 1972, assuming the satellite particle distribution to be described 
by the generally accepted relationship = 2,5R^. Figure 4 summarizes 
the cross section measurements for various ions in the energy range 1 keV 
to 200 keV in the form of normalized charge exchange lifetimes Thus., 

using these figures the equatorial lifetimes can be determined for any of 
these ions at a specific energy and L-value. The charge exchange lifetime 
for the mirroring particles with a specific mirror latitude can be deter- 
mined by t = t cos " or by using equation (7). 

The results of this paper are based on our current knowledge of the 
planetary atmosphere neutral hydrogen distribution and the charge exchange 
cross sections for various ions. The charge exchange lifetimes of the 
equatorial particles are much shorter and the dependence on mirror latitude 
is much weaker than the Liemohn (1961) results. Thus, we recommend that 
the ionic composition of the ring current during the recovery phase of 


22 



magnetic storms should be re-examined in light of these current results. 

To highlight this conclusion, we show, as an example, how these 

results affect the particle flux decay as a function of mirror latitude. 

Assuming at time t=o, an isotropic and unit flux distribution for 10 keV 

0"^ and H’’” ions at L = 4, and assuming geophysical parameters corresponding 

to the February 24, 1972 magnetic storm, we show in Figure 5 the mirror 

latitude distributions as predicted by the new lifetime values after time 

periods which are multiples of 7.8 hours (roughly the orbital period for 

Explorer 45). The dotted lines show the distribution predicted by using 

Liemohn's mirror latitude dependence of t =t cos^X and the solid lines 

mem 

3 5 

correspond to the dependence of cos " Xj^i for both cases the equa- 

torial lifetimes, are the new values summarized in this paper. 

It is clear that the flux distributions as a function of mirror 
latitude decay much slower using the new functional dependence.^ Two 
additional points are also evident from this figure. First, assuming 
comparable initial fluxes, the 0 ions after about one day would exhibit 
a pitch angle distribution significantly separated from the H ions which 
would have more rapidly decayed. This can aid in the determination of 
ion composition in the ring current. Second, the loss cone for these 
particles is about 60° mirror latitude and, therefore, these distributions 
will be quite flat with sharp cut-offs near the loss cone for many hours 
after on initial isotropic distribution. Obviously the ions with longer 
equatorial lifetimes (eq. 0^ and He'*’) will exhibit this shape for a 
longer time. 

With these new lifetime values it will be Indeed a challenge to 
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re-examine some of the previous results, and exciting to interpret the 
new measurements and to gain a more complete understanding of our ring 
current environment. 
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FIGURE CAPTIONS 


Figure 1.. 


Figure 2: 


Figure 3: 


Figure 4: 


Figure 5: 


Functional dependence of atomic hydrogen density on the 
radial distance from the center of earth as given by Johnson 
and Fish model (shown as dashed line) and the Chamberlain 
model (shown as solid line). A unit density and tempera- 
ture of 1200°K at the exobase Tc are assumed. For Chamberlain 
model, values of = 2.5 is used. 

Functional dependence of atomic hydrogen density distribu- 
tion as given by the Chamberlain model for a range of 
exobase temperature. Panels (a), (b), (c) and (d) correspond 
to = 900°K, 1000°K, 1100°K and 1200°K respectively. For 
each panel, dashed line corresponds to all satellite par- 
ticles included, solid line corresponds to partial 
population of satellite particles (R = 2.5 R ), and dash- 
dotted line corresponds to no satellite present. Unit den- 
sity at the exobase is assumed. 

Change in neutral hydrogen density as a function of radial 

distance (in Earth radii). Variation in exobase temperature, 

T^, of 900°K to 1200^K, and in density of 4 x 10^cm~^ to 1.1 
5 -3 

X 10 cm is assumed. 

Normalized charge exchange lifetimes (1/av) in seconds per 
cm for various ions in the incident energy range of 1 to 200 keV. 
Evolution in pitch angle distribution as a function of mirror 
latitude for charge exchange decay, starting with a unit 
isotropic distribution at time t 0. Each panel corresponds 
to a time interval of 7.8 hours. Solid lines show the decay 
as predicted by t /t = cos dependence and dotted line 

iM 6 m 

shows the decay as predicted by 'c^/Tg= dependence. 
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